Causality and the Doppler Peaks 
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A considerable experimental effort is underway to detect the 'Doppler peaks' in the angular power 
spectrum of the cosmic microwave anisotropy. These peaks offer unique information about structure 
formation in the universe. One key issue is whether structure could have formed by the action of 
causal physics within the standard hot big bang, or whether a prior period of inflation was required. 
Recently there has been some discussion of whether causal sources could reproduce the pattern of 
Doppler peaks produced by the standard adiabatic theory. This paper gives a rigorous definition 
of causality, and a causal decomposition of a general source. I present an example of a very simple 
causal source which mimics the standard adiabatic theory, accurately reproducing the behaviour of 
the local intrinsic temperature perturbations. 



Existing theories of cosmic structure formation are of 
two types. In the first, the hot big bang is assumed to 
have started out smooth. Structure then forms as the re- 
sult of a symmetry breaking phase transition and phase 
ordering. In the second, an epoch of inflation prior to 
the hot big bang is invoked. Whilst both mechanisms 
are causal, causality imposes a much stronger constraint 
in the former case (Figure 1), because the initial con- 
ditions for the perturbation variables are established on 
a Cauchy surface £ within the hot big bang (Figure 1). 
The causal nature of the Einstein-matter field equations 
then implies the vanishing of all correlations between all 
local perturbation variables at spacetime points whose 
backward light cones fail to intersect on E. In the infla- 
tionary case, by construction the relevant surface £ lies 
so far before r = that there is no useful constraint. 
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FIG. 1. Causal constraint on theories of structure forma- 
tion where the standard hot big bang starts out homogeneous. 
The vertical axis shows conformal time r, with r = corre- 
sponding to the initial singularity (in the absence of inflation) . 
Correlations between any local variables at any two spacetime 
points vanish if their backward light cones fail to intersect on 
E, the spacelike hypersurface r = tpt just prior to the phase 
transition. In inflationary theories, there is no singularity at 
T = 0, instead there is a preceding epoch of inflation during 
which longer range correlations are established. 

Could observations distinguish the causally con- 
strained theories from inflationary ones? The cosmic mi- 
crowave anisotropy is the best hope of a direct probe, 



giving us a picture of the universe on the surface of last 
scattering. This surface cuts through many regions which 
were 'causally disconnected' (quotes indicate a standard 
hot big bang definition) at that time. If the only con- 
tributions to the microwave anisotropy were local ef- 
fects, like temperature and velocity perturbations in the 
photon-baryon fluid, one could check whether 'superhori- 
zon' perturbations were present by measuring the auto- 
correlation function of the anisotropy map. If this was 
consistent with zero beyond some angular scale (twice 
that subtended by the 'causal horizon' at last scattering, 
of order 2° with standard recombination), one could con- 
clude that the perturbations were indeed causally con- 
strained. 

The complication that spoils this test is that a signif- 
icant component of the microwave anisotropy is gener- 
ated after last scattering, by the integrated effect of time 
dependent gravitational potentials along photon paths. 
This is after all how cosmic defects produce a scale in- 
variant spectrum of microwave anisotropies on very large 
angular scales (consistent with the COBE results) even 
though these theories are causally constrained. 

Nevertheless, the local contributions to the microwave 
anisotropy do have a signature distinguishing them from 
the foreground due to the integrated effect. This is 
the presence of 'Doppler' peaks in the angular power 
spectrum, caused by phase-coherent oscillations in the 
photon-baryon fluid prior to recombination. The loca- 
tion of these peaks is mainly determined by the temper- 
ature perturbations in the photon-baryon fluid, a com- 
pletely local effect. This Letter will address the question 
of whether the peak locations can be used as a discrimi- 
nator between inflationary and non-inflationary theories 
of structure formation. 

Crittenden and I suggested a connection between 
causality and peak location following an analysis of 
the cosmic texture theory, in which the Doppler peaks 
are phase-shifted relative to those in standard inflation, 
the biggest peak occurring at higher multipole I (smaller 
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angular scales) than in the standard inflationary theory. 
Albrecht, Magueijo and collaborators (2| raised the im- 
portant issue of decoherence, and gave a detailed discus- 
sion of the behaviour of the Doppler peaks for different 
models of causal sources, in particular those motivated by 
the study of Robinson and Wandelt ||] . Most recently Hu 
and White Q have addressed these issues, claiming that 
the pattern of Doppler peaks predicted by the simplest 
inflation model is 'essentially unique and its confirmation 
would have deep implications for the causal structure of 
the early universe'. In this letter I develop a formalism 
for dealing with decoherent but causal sources. I ex- 
hibit a causal source which closely mimics the inflation- 
ary pattern of temperature perturbations in the photon- 
baryon fluid, and the corresponding contribution to the 
microwave anisotropy. 

Structure formation within the standard hot big bang 
requires the presence of a source term in the Einstein 
equations in addition to the usual metric and matter 
variables. Cosmic string and texture each provide an 
example of such a sources. The perturbations are most 
simply dealt with in the fluid approximation, which is 
reasonable for our purposes. In the synchronous gauge, 
the relevant equations read |J 

5 C + -Sc = 4^^(1+34^ + 5, (1) 

N 

S R + -(1 - 34)4 = 4v 2 s R + U'c + ~(1 - 34)<5 C 
a o o a 

(2) 

Dots denote derivatives with respect to r, a(r) is the scale 
factor, 8c and Sr are the contrast in dark matter and ra- 
diation densities, and cs is the speed of sound. The sum 
over N includes dark matter, the photon-baryon fluid 
and neutrinos. I assume the 'canonical' parameter values 
Q.CDM = 0.95, VL B = 0.05, and h = 0.5. The fluctuating 
part of the external source is taken to have stress energy 
tensor 0^, and S — 47rG(9oo + ©«)• The initial condi- 
tions to be used with ([!]) and (|J) are the vanishing of the 
pseudoenergy too = ©oo + Dw Pncl 2 5n + {a/ a)5 c / {^G), 
and adiabaticity, Sr = S v = f^c- This corresponds to 
starting with a universe in which the energy density and 
space curvature are uniform. 

It is a good approximation to treat the radiation 
as a perfect fluid, obeying the equations ([!]) and (||) 
from initial conditions set up deep in the radiation 
era, up to recombination. The intrinsic temperature 
perturbation from a Fourier mode k is then given by 
(5T/T)j(k) = j5R(k), and its contribution to the an- 
gular power spectrum of anisotropies is given by C; cx 
J k 2 dk(ST/T) 2 (k)ji(kT ) 2 with r the conformal time to- 
day This dominates the anisotropy pattern on small 
angular scales. 



Equations (|l|) and (g) are linear, and it follows that 
all correlations between local observables are completely 
determined by the unequal time correlation function of 
the source stress energy tensor. In particular, for S the 
causality constraint reads 

£(r,r,r') =< S(r,T)S(0,T') >= r > r + r' (3) 

The sharp edge on £ leads to oscillations in its three di- 
mensional Fourier transform £(fc, r, r') at large k. Inte- 
gration by parts produces 

R R ' 

£{k, t, t') ~ [--^cosfcr + j± sinfcr + ..\l +T , (4) 

where R(r) = r£(r), and R. r = dR/dr. (If £ ~ r~ 2 at 
small r, as it does for strings, then £ has an additional fc -1 
term) . The leading term is not necessarily oscillatory, but 
there must be oscillatory subleading terms. Most of the 
ansatzes for £ in the literature do not have this feature 
and are therefore manifestly acausal. They may still be 
useful as approximations, but it is desirable to develop a 
formalism in which causality is rigorously built in. 

The discussion simplifies if we assume scaling |5| . Then 
dimensional analysis implies that 

£=< S(k,T)S*(k,T') >=T-h , -ix(kT,kT'). (5) 

Regarded as a matrix with indices r, and r', X is real 
and symmetric and can therefore be represented as: 

X(kT,kr') =^P a / a (fcr)/ a (fcr / ) (6) 

a. 

with f a {kr) a set of orthonormalised eigenfunctions of 
X(kT,kr') regarded as an integral operator, with cor- 
responding eigenvalues P a . In the terminology of [jy, 
fa(kr) is a 'master' function. As in quantum mechanics, 
we have a pure, 'coherent' state if P a is nonzero for only 
a single value of a, otherwise we have a mixed, 'incoher- 
ent' state. Equation (||) shows that a general source may 
be represented as an incoherent sum of coherent sources. 
The P a 's must be positive for all a because they are the 
expectation value of a quantity squared. 

This representation is useful because the contribu- 
tion of each individual term in the a sum is straight- 
forwardly calculable, by using the source t~? f a (kr) in 
the linearised Einstein equations. A bonus is that the 
assumption of scaling allows one to infer the correct 
initial conditions for the perturbations. For small fc, 
@ assures us that, if £~(0) exists, f(kr) must tend to 
a (possibly zero) constant. Then energy conservation, 
6oo + (a/a) (0oo + 0«) ~ for fcr << 1, and the as- 
sumption of scaling (by dimensions, 0oo « r^) allow 
one to unambiguously determine the contribution to 0oo 
appropriate to each f a (kr). 

Now let us return to the causality constraint (^) . Con- 
sider a single term in the sum over a in (^) . The contri- 
bution it makes to < S(r, t)S(0, t') > is proportional to 
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the convolution of f a (r, t) with f a (r, r'). If the f a (r,r) 
have compact support, a simple argument Q shows that 
fa(r,T) = for all r > r. This gives a nice geometrical 
picture of how the causality constraint works. For each 
a the master function f a (r, r) is the profile of a ball of 
radius r, and the convolution of f a {r, t) with f a (r,T') 
clearly vanishes if the separation of the ball centers is 
greater than r + r'. 

Determining the form of the f a and P a would be very 
interesting in any particular causal scenario. Here how- 
ever, I want to see whether anything useful can be learnt 
by considering all possible / Q 's and P a 's. The power 
spectrum in the general case is is just a sum of the power 
spectra for different such /'s with positive coefficients, so 
if for example we can show that the C; for every / has 
positive slope for I < l m ax, it follows that the total power 
spectrum will too. In this way we can set a lower limit 
on the location of the first Doppler peak. 

A basis for all functions /(r) is provided by the fam- 
ily r2 f( r , T ) = <>( r ^ At), with < A < 1. In Fourier 
space we have f(k,r) = sin Akr / (Akr) . In one extreme, 
with only short range correlations, f(k,r) is nearly con- 
stant. In the other, it has its first zero at kr = n. The 
equal time correlation functions corresponding to this 
family of master functions are smooth functions: one find 
£(r,r,T) - l/(rr 3 ) for r < 2At, £(r) = for r > 2At. If 
the master functions f(r) do not change sign for all r < r 
(which is unlikely), then any f(r) can be represented as 
a sum of the above basis functions with positive coeffi- 
cients. 

We now proceed to solve equations ([!]) and for this 
family of source functions, with S = /fe(r)/r2. Each 
Fourier mode of 5r starts out small and grows (like rs). 
After horizon crossing it oscillates as an acoustic wave. 
At the 'instant' of last scattering, all modes are caught at 
a particular phase of their oscillations, and those which 
are at maximum amplitude produce the Doppler peaks in 
Ci . Figure 2 shows the time evolution of a single Fourier 
mode of 5r and 5c in the two extreme cases (A = 0, 
denoted O, and A = 1, denoted X), and in the stan- 
dard adiabatic theory with no source. The approximate 
scale invariance means that the same graph very roughly 
represents 5R{k,T rec ) as a function of k at recombina- 
tion. One can translate fcr rec into multipole moment I 
by the approximate relation I ~ kr® ~ b0kr rec . Peaks 
in 5r(A:,t) 2 are, through the integral given above, trans- 
lated into peaks in G\. 

In the causal theories, 5 c is forced to start out grow- 
ing with sign opposite to the source iS, because the total 
pseudoenergy Too must initially be zero. As time goes on, 
S starts to drive 5c- 

US always has the same sign, as in the case A << 1,5c 
changes sign as it becomes driven by S. The forcing term 
for the radiation, 5c, then changes sign around kr = 1, 
so while 5r initially grows with the opposite sign to S, 
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FIG. 2. The evolution of perturbation modes for radiation 
5r and dark matter 5c as a function of conformal time r, for 
fixed k. The value of k chosen has kr rec = 5, and in the infla- 
tionary theory the corresponding mode has its first maximum 
at recombination. The superscript A denotes the standard 
inflationary theory, X shows the extreme causal model which 
mimics the case A, and O shows the opposite extreme causal 
model which mimics instead the texture model. 

it is later driven to the same sign as S. Because the sign 
change in 5c occurs early (at fer ~ 1) the first oscillation 
in the radiation has small amplitude. Because Ci is really 
an integral over fc, as mentioned above, the first Doppler 
peak, at I ~ 120, is smeared by the contribution of higher 
k, and may be effectively 'hidden'. The main peak is that 
due to the next oscillation, the one that is really 'driven' 
by S. This one occurs at / ~ 380, compared to the main 
peak in the standard adiabatic case at I ~ 220. Inside the 
horizon, the radiation oscillates sinusoidally, and higher 
peaks occur at shifts Al « 280n, n — 1,2,3,.. to the 
right. It is interesting that this case (A << 1) reproduces 
the main features of the texture models presented in Q| 
and§. 

Next, consider the case where there is a sign change 
in >S around horizon crossing. As before 5c and 5r start 
out with the opposite sign to S, because of compensation. 
But here, if S changes sign early enough, 5c does not have 
to change sign. The radiation forcing term 5c is always 
positive, and the first peak in 5r is not small. As can 
be seen in Figure 2, the extreme case A = 1 mimics the 
standard adiabatic model rather closely. I have computed 
the power spectrum 5R(k,r) 2 at recombination in the 
A = 1 theory, for Q, B = 0.05, 0.1, and 0.2 (and VLcdm = 
1 — VLb)- In all cases the result is similar to the analogous 
standard adiabatic theory, both in the peak location and 
in the pattern of peak heights. 

It may be useful to visualise this in terms of the C\ 
spectra. Figure 3 contrasts the standard adiabatic C;'s 
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with those of the texture model 0,0. The simple family 
of causal theories I have just discussed roughly speaking 
interpolates between these two curves. At A = 1 the first 
peak is close to that of the adiabatic theory, and as one 
decreases A, the peaks move to lower I, The first peak 
decreases in amplitude, for the reason discussed above, 
and moves leftward from I ~ 220 to I ~ 120. The second 
increases in amplitude, and moves from I ~ 500 to I ~ 
380. 




1000 

FIG. 3. The anisotropy power spectra for the standard in- 
flationary theory (dashed line) and the texture theory (solid 
line) . The family of sources studied here produces peak loca- 
tions which approximately (at the ten per cent level) interpo- 
late between these two cases. As the parameter A is dialed 
from 1 to 0, the position of the first peak moves from I ~ 240 
to I ~ 120, and decreases in amplitude. The higher peaks 
shift down in I by a similar amount, with the second peak 
growing to become the highest peak. Note that the texture 
curve shown here includes the vector and tensor contributions, 
which help to emphasise the 'hidden' peak at I ~ 120. 

The calculations reported here include only the Sr con- 
tribution to G\ , whereas the theoretical curves shown in- 
clude all contributions, including the radiation velocity 
terms and the Sachs- Wolfe integral. In the inflationary 
case, the Sr term alone gives a good overall represen- 
taion of the Ci spectrum for I > 100, and particularly 
the peak locations. It is therefore reasonable to expect 
that, at least in a similar situation where the Sachs- Wolfe 
integral and radiation velocity terms are sub-dominant, 
the complete spectrum of Cj 's will be similar. Computa- 
tion of the Sachs- Wolfe integral, however, requires spec- 
ifying the full stress tensor 6^„ of the source, not just 
the combination S which was sufficient here. In a com- 
panion paper [jll[, I extend the definition of the source 
stress tensor to one in which the Sachs- Wolfe integral is 
small, and show that with this extended definition of the 
source, the qualitative conclusions reached here hold for 
the complete C\ spectrum. 



Continuing with the discussion of the Sr contribution 
alone, what general conclusions can be drawn? Any mas- 
ter function / can be represented by a sum of the ba- 
sis above: since 6r follows the same t' evolution up to 
kr ~ 2.5 for all of them, it follows that this will be true 
in the general case. Translated into C[s, this means that 
they cannot have a peak below I ~ 120. This is the 
real, and perhaps disappointing, causality constraint on 
scaling causal sources. Can we push the first peak to 
higher I than in the standard theory? Within the family 
considered, the limit for the first peak is close to the adi- 
abatic position, I ~ 240. It follows this is an upper limit 
any theory where the master functions f(r/r) are strictly 
non-negative. However, if the f(r/r) do change sign, the 
first peak may be pushed to much higher I. For example 
adding the negative of the 'O' case to the 'X' case in Fig- 
ure 2 pushes the first Doppler peak to I ~ 400. Similar 
examples produce first Doppler peaks at even higher I. 

In conclusion, I have proposed a new formalism within 
which causal sources can be studied. As a first ap- 
plication, I have exhibited a simple family of strictly 
causal sources which produce local contributions to the 
anisotropy which approximately interpolate between the 
standard adiabatic prediction and the texture prediction 
given in jlj and . In a companion paper JTT| the source 
model is extended to define the full stress tensor 0^,,, 
in such a way that the Sachs Wolfe integral is a sub- 
dominant contribution to the anisotropy. This extended 
model shows that the causality constraint cannot on its 
own be used to distinguish between inflationary and non- 
inflationary theories, and extra details are needed. It will 
be very interesting to determine the master functions f a 
and coefficients P a for specific scenarios. 
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PHY90-21984, and the David and Lucile Packard Foun- 
dation. 
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